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SUMMARY 


Curved pipes have unique stress distributions and flexibility 
properties which cannot be easily reproduced using the beam elements 
of the NASTRAN finite element structural analysis program. The purpose 
of this research effort was to enhance the application of NASTRAN to 
the structural analysis of piping systems by introducing a curved pipe 
element capability. The theory used for the element stiffness matrix 
was developed from an existing flexibility matrix. The curved pipe 
capability was implemented using the NASTRAN user dummy element. 

The dummy element approach is compatible with most of the NASTRAN 
rigid formats and requires less sophisticated programming than needed 
for incorporating a new element into NASTRAN. The similarities between 
the NASTRAN bar element and the curved pipe element made it possible 
to use the bar element subroutines as the basis for the necessary 
curved pipe subroutines. An existing program was used to link the 
subroutines to NASTRAN. Test problems have been prepared and input 
and output pertaining to the curved pipe element are described along 
with solution results. 
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INTRODUCTION 


Object and Scope 

Curved pipes have unique structural properties under flexure, which 
cannot be easily modelled using standard NASTRAN elements. Flexure causes 
an ovalization of the cross section, resulting in a significant change 
in the stress distribution in the curved pipe. Longitudinal bending 
stresses in the extreme fibers are relieved and high circumferential 
stresses occur. As a result, the peak stress is in the circumferential 
direction. This peak stress is conventionally related to that which would 
exist in a straight pipe subjected to the same moment through a "stress 
intensification factor." An additional effect of the altered stress 
redistribution is a decrease in the stiffness of the curved pipe. The 
ratio of the stiffness of a straight pipe to the stiffness for the curved 
pipe is denoted the "flexibility factor" and is related to the geometry 
of the pipe. This behavior could be modelled using NASTRAN bar elements 
by: (1) using a sufficient number of bars to model the pipe curvature; 

(2) accounting for stiffness reduction by reducing the moments of inertia; 
and, (3) using appropriate stress recovery factors to simulate the stress 
intensification factor. This approach would be cumbersome for a piping 
system with more than a few bends. The object of this research effort is 
to simplify the analysis of piping systems using NASTRAN by developing 
a curved pipe element. This element would minimize the effort of the NASTRAN 
user and would maintain as many of the standard NASTRAN capabilities as 
possible. The NASTRAN user dummy element approach was chosen for 
implementing this new capability because it is the simplest approach to 



developing a new element and because the present restricted capabilities 
of the curved pipe element definition make it unsuitable for permanent 
incorporation into NASTRAN at this point. Also, permanent incorporation 
should be preceded by a thorough checkout of the proposed element using 
the dummy element approach. The dummy element can be used for all of the 
NASTRAN structural rigid formats except the piecewise linear analysis. 

The element is implemented by first writing FORTRAN IV subroutines for 
the portions of NASTRAN which are dependent on the element used and then 
linking these subroutines into NASTRAN. Capabilities which account for 
offsets of bar ends from grid points, released constraints at bar ends, 
temperature loads and enforced deformation loads are included for the 
element. No differential stiffness matrix is generated. A lumped mass 
approach is used to described inertia properties. Input and output are 
similar to that for standard NASTRAN elements. 

A considerable amount of computer programming was required in 
order to develop the necessary subroutines for the curved pipe element. 
Maria V. Stephens of the Engineering Analysis Branch, NASA Langley 
Research Center, provided this service and also provided much needed 
assistance in dealing with the complicated programming logic involved 
in the general NASTRAN program. Her assistance is gratefully acknowledged. 

Also, this research effort was conducted as part of the NASA-ASEE Summer 
Research Program. 
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Literature Review 


Although several finite element programs include a curved pipe 
element, details about how the elements function are unavailable. The 
basis for any finite element is a stiffness or flexibility matrix which 
relates the forces acting on the element to the resulting displacement. 
Although several sources (References (1), (2), (3), (4)) have developed 
and used stiffness matrices for curved beam elements, which have 
stiffness properties similar to those of curved pipes, relatively little 
has been done specifically for curved pipes. A flexibility matrix 
developed in Reference (5) and modified in Reference (6) is used as 
the basis for this analysis. 
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ANALYSIS 


Development of Subroutine Relationships 

In order to introduce a new element capability to NASTRAN via the 
dummy element approach, five FORTRAN IV subroutines must be developed 
which: (1) compute the stiffness contributions of the element to the 

adjacent grid points; (2) compute mass contributions to the adjacent 
grid points; (3) compute thermal and enforced deformation loads; 

(4) set up force recovery matrices; and, (5) recover forces and stresses 
for output. These subroutines must then be linked to the NASTRAN program 
before the dummy element can be used. The theory involved in these 
subroutines is discussed in this section. 

The development of the subroutine relationships is complicated 
somewhat by the different coordinate systems involved. NASTRAN allows 
the user to define displacements at each grid point in a coordinate 
system of the user's choosing. This set of coordinate systems is then 
used by NASTRAN in solving for the grid point displacements. However, 
the properties of an element are generally best defined using a coordinate 
system peculiar to that specific element. These two coordinate system 
sets are related to each other in that both coordinate system sets are 
defined in terms of a single coordinate system, called the "basic" 
coordinate system. The relationship between the displacement coordinate 
system and the basic coordinate system is defined on bulk data cards; 
information needed to define the element coordinate system in terms of 
the basic system is inputted on bulk data cards, but the actual relation- 
ships are developed in the first of the five "dummy" subroutines. 
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The theory used in developing the stiffness matrix for the curved 
pipe element was based on the flexibility matrix described by Chen 
Reference (1). The geometry and notation used by Chen is shown in 
Figure (1). (The "z" direction is out of the paper). Chen s matrix 
is written with reference to the center of the arc described by the 
pipe centerline. It is assumed that end (b) of the bar shown in 
Figure (1) is fixed and that a rigid connection is made between end 
(a) and the center point 0, so that the flexibility matrix relates the 
deflections at point 0 to forces applied at point 0, or, mathematically, 


where 


<u_> ' fu 0 U 0 U 0 ®0 9 0 9 0 } 1S a 6 * 1 C ° 1U,nn "' atr1X 

0 u y u y Z X y z 


giving the 3 component deflections and 3 component rotations at point 

0 with reference to the coordinate system xyz. The flexibility matrix 

[f] is a 6 x 6 matrix of influence coefficients. The force vector, 

,r , _ rr c C M M M } is a 6 x 1 column matrix giving the 3 
4 ° } ' °x °y °z °x °y °z 

component forces and the 3 component moments at point 0 with reference 
to the coordinate system xyz. Equilibrium conditions for the element, 
the reciprocal theorem, and Equation (1) can be used to get relationships 
between forces and displacement at the two ends of the element. Specifically, 
the following relationships are required: (1) a stiffness matrix K ag , 

relating displacements at end (a) to forces at end (a); (2) a stiffness 
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Figure (1). - Geometry and Notation for Chen Flexibility Matrix 
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matrix, K^, relating displacements at end (a) to forces at end (b); 

(3) a stiffness matrix, K ^ relating displacements at end (b) to forces 
at end (a); and (4) a stiffness matrix, K bb , relating displacements at 
end (b) to forces at end (b). The notation used for these forces and 
displacements is given in Figure (2), with the z direction (out of the 
paper) omitted. 

It is convenient to work with Equation (1) in its inverted form. 


{F o } 


6 x 1 


[K°] {u n > 

6 x 6 0 6 x 1 


where 

[K°] = [f]" 1 . 


( 2 ) 


Since (o) - (a) is a rigid connection in Figure (1), forces applied at point 
(o) could be transferred to end (a) of the element by the equation 


where 


{F a } 


6 x 1 



] 

6x6 


{F o } 


(3) 


[B 


oa 


] = 


I 

0 0 0 
0 0 R 

0 -R 0 



Displacements at point (0) can be written in terms of the displacements 
at end (a) using the equation 


< u „> . - t B oa>; . fu a } 

6x1 6xo 


(4) 
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Figure (2). - In Plane Forces and Displacements 
for the Curved Pipe Element. 
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Equations (2), (3), and (4) can be combined to get a relationship 
between forces applied at end (a) and displacements at end (a): 


(F> = [K aa ] 

a 6 x 1 aa 6 x 6 


{u a } 


6x1 


( 5 ) 


where 


[K„] - IB ] [K°] 

aa oa 6 x 6 6x6 


[B ] 
L oa J 


( 6 ) 


Equilibrium requires that forces applied at end (a) be balanced 
by equal and opposite forces at end (b), in the absence of intermediate 
loads on the element. Thus, the reaction forces at end (b) associated 
with forces applied at end (a) can be found by transferring the forces 


at end (a) to end (b), redefining them in terms of the coordinate system 
at end (b), and changing the signs. This process can be written in 


matrix notation as 


where 


<F h ) ’ - !B a J (FJ 

b 6 x 1 >6 


t B ab ] * 


cos 4* 
-sin 
0 
0 
0 

R sin 


sin \p 
cos ip 

0 

0 

0 

R( 1 - cos i|j) 


6x1 


0 0 0 

0 0 0 

1 0 0 

-R sin cos \p sin ip 

R( 1 - cos \p) -sin cos ij; 

0 0 0 


(7) 

o’ 

0 
0 
0 
0 
1 
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Equations (5) and (7) can be combined to yield 


(F b> , = l K bal '“a 5 , . 

6x1 6x6 6x1 


where 


tK ^ ] 6 x 6 = " ^6 x 6 [K - ] 6 x 6 


( 8 ) 

(9) 


Forces at both ends of the element due to displacements at end b 
could be found using the same approach as for displacements at end a. 
This is not necessary when a relationship between the displacements at 
end (b) and the forces at end (a) are desired, since the reciprocal 
theorem requires that 


{F a> 


6 x 1 


[I W 


6x6 


<%) 


6 x 1 


( 10 ) 


If end (a) is assumed fixed while end (b) is assumed rigidly attached 
to point (0), displacements at end (b) can be related to the forces 
applied at end (b) (using the coordinate system indicated in Figure (2)) 
by the relation 


lF t> , - ' K bbi s <V 

6x1 6x6 


(ID 


where 


t^b 1 - l B 0 b ] [K ° ] [B ob ,T 

DD 6 x 6 0D 6 x 6 6x6 00 6x6 
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and 


B .] 

ob 


cos ip 
-sin \p 

0 

0 

0 

R sin ip 


sin rp 
cos ip 

0 

0 

0 

-R cos 


0 

0 

1 

0 

R 

0 


0 0 0 

0 0 0 

0 0 0 

cos (Jj sin ip 0 

-sin \p cos ip 0 

0 0 1 


The relationships between forces and displacements at the ends of the 


element can now be summarized by combining Equations (5), (8), (10), 
and (11) to get 


where 



= [K E ] 


12 x 12 



<12 x 1 


[K E ] 


12 x 12 





I 


T 

ba 


bb 


J 12 x 12 


( 12 ) 


( 13 ) 


At times it is desirable to disconnect various degrees of freedom of the 
ends of the element from the adjacent grid pont. This capability is 
included in the curved pipe element through the NASTRAN pin flag routine. 
This routine accounts for disconnected degrees of freedom by operating 
on the element stiffness matrix, [K^]^ x ^ to P rocluce "reduced" 

p 

stiffness matrix [K ] ^ x ^2' 

At this point forces at the ends of the element have been related 
to displacements at the ends of the element, as defined by the coordinate 
systems shown in Figure (2), by the equation 


11 



{F} 12 x 1 = [I(R] I2 x 12 {u} 12 x 1 


where 


[K R ] 


12 x 12 


ba bb 

The equation solved in NASTRAN is of the form 


„R v R 

a a Da 


K 


R 


,R 


(S> = [K G ] (u G ) 


(14) 


(15) 


where 

{S } includes all external loads applied to the grid point, 
and is therefore equal to the sum of the internal forces 
applied by the elements to grid point, 

p 

[K ] is the global stiffness matrix, and can be partitioned as 



(u G > is the complete set of all grid point displacements. 

p 

Both matrices (S } and { u } are defined in the displacement coordinate 

p 

system at each grid point. The global stiffness matrix [K ] is an 
assemblage of the stiffness contributions of all the elements. Thus, if 
the force at the adjacent grid point can be written in terms of the force 
at the end of the element and the displacement at the end of the element 
can be written in terms of the displacement at the adjacent grid point, 
the element's contribution to the global stiffness matrix can be determined 
using Equation (14). The coordinate systems involved in determining this 

12 
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contribution is shown in Figure (3). Offset vectors (tu^ = 

h z * “A 3 \ and = “Bj \ + “B 2 *b 2 + "b 3 »b 3 > are also shown - The 

displacements at end (a) of the element can be written in terms of the 
deflection at grid point A using the relation 


J a^ . ~ ^eb r 

6x1 6x6 


[C a ][ E a ] (u a ) 

M A 6 x 6 H 6 x 1 


where 



and determines the displacement at (a) in the displacement 
coordinate system at grid point (A) given displacements at (A); 

[C.] is an orthogonal transformation matrix which redefines 

M 

in the basic coordinate system, vectors defined in the displacement 
coordinate system at A; 

(a) ..... 

[T . ] is a transformation matrix which redefines in the basic 

1 eb 

coordinate system, vectors defined in the element coordinate system 
at end (a); since this is an orthogonal matrix, the inverse (which 
redefines in the element coordinate system at end (a), vectors 
defined in the basic coordinate system) is equal to the transpose. 
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Figure (3). - Coordinate Systems Used in Element Subroutines. 
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A force vector at grid point A can be written in terms of an equivalent 
force vector at element end (a) using the relationship 


<V = [E/ [C/ [T^HF.) 


( 17 ) 


Combination of Equations (17), (5), and (16) yield 


{F A } = [K AA ] {u A } 

A 6x1 ^6x 6 o x 1 


where 


(a) 


=(l T eb'> ' C A> 1 E a 1> T f K a a ] <t T eb l t C A» E A]> 


(a) 


6x6 

Similar efforts produce the matrices 


(18) 


(b) 


UT e '] [C B ] [E b ]) T [K^ a ] ([T eb ] [C A ] [E A ]) 


(a) 


G 

[K BA ] C c ‘ Ul eb J LU B J Ll -B J/ LIN ba J VL 'eb 

0 X u 


(19) 


and 


(b) 


(b) 


[K BB ] c = ([T eb' ] f £ B^ T [K bb^ [T eb ] t C B ] [ £ B ] ) 

6x6 


(20) 


where 


W - 


I "lllr 


■d)r 


"0)r 


0) r 


6x6 


(b) 

and [C D ] and [T . 3 are transformation matrices at B equivalent to matrices 
B 6D 

(a) 

[C A ] and [T eb J at A. 
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Finally, application of the reciprocal theorem yields 


(k ab ] 


6x6 



T 

] 

6x6 


( 21 ) 


The submatrices given by Equations (18), (19), (20), and (21), provide 
the contribution of the curved pipe element to the total global stiffness 
matrix of the structure. 

A lumped mass approach was used to define the inertia properties 
of the curved pipe element. It is assumed that the total mass of the 
element can be divided equally between the element ends. Second mass 
moments are not included in the mass matrix. This omission should have a 
negligible influence for curved pipes with small arc but could be significant 
for curved pipes with larger arcs. This is because of the lack of 
compensation for the offset of the center of gravity from a tangent 
connecting the ends of the element. If it is desired to include the 
effects of the second mass moments, NASTRAN has a capability for directly 
inputting these properties. Forces applied to either end of the element 
can be related to the corresponding acceleration components by the equation 


(F } = [Ml (u> 

a 6 x 1 a 6 x 6 a 6 x 1 


where 


{u a ) is a 6 x 1 column matrix giving the linear and angular components 
of acceleration of end (a) of the element. 
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and 


[Ml 

a 


m 

2 


m 

2 


= f M b ] 


( 22 ) 


m 

2 


I — 


with 


m 


= R * {p J (d 2 - (d - 2t) 2 ) + y} 


p = mass density of pipe material 
y = non-structural mass per unit length 

The NASTRAN solution process requires a relationship between the forces 
and accelerations at the grid points rather than the ends of the element. 

A process similar to that used for the stiffness matrix yields 


where 


{F a } = [m aa ] {u a } 

A 6 x 1 M 6 x 6 A 6 x 1 


(a) 


(a) 


‘ M A«l c ■ <' T eb 1 'V f E A ! > '"a 1 < [T eb > 'V ' E A» 

D X 6 


Simi larly , 


(F r > = [Mor] {u fl } 

B 6 x 1 BB 6 x 6 A 6 x 1 


where 


f (b) T T (b) T 

[N B b 1 = ([T eb) rc e l t E B ]) T [MbJ (t T eb ^ IVM 

OXO 


( 23 ) 


( 24 ) 
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This formulation assumes no direct relationship between forces at 
(A) and accelerations at (B) and vice versa. Thus Equations (23) and (24) 
define the total mass contribution of the element to the structural system. 

Stresses due to thermal expansion loads are significant in piping 
system analysis. Stresses due to misfit may also be of interest. A 
capability for handling both types of problems has been included in 
the curved pipe element theory. It is assumed in both cases that the 
ends of the element move together or apart, whichever the case, without 
any rotation occurring. The forces required for this compression or 
extension are calculated and added to the existing external loadings at 
the grid points. Disconnected degrees of freedom are taken into account 
in determining these forces. Once the problem has been solved, the forces 
which would have resulted from the enforced deformation are calculated and 
subtracted from the forces calculated from the movement cf the ends of the 
element. 

The element is assumed to be fixed at end (b) (see Figure (3)) 
so that the movement associated with the applied deformation or temperature 
expansion can be interpreted as a displacement of end (a) in the plane of 
the element. The equivalent grid point loads are 

lF A> , . ■ tc/ [Teb’l [K a R a ] (0L) (25a) 

oxl 6x1 

and 

< F B> , , - [E B ,T (C B ]T I T eb ] [<,) fDL a ) (25b) 

oxl 6x1 
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where 


{DL } = the displacement at end (a) in the coordinate system 
d 

at end (a). This displacement is automatically calculated in the 
subroutine given the enforced deformation and temperature information. 

The element temperature is calculated from the relation 

T =l< T o + 2T m + T d 

where 

T = the temperature used to calculate expansion 
T q = the temperature at the outside of the pipe wall 

T = the temperature at the middle of the pipe wall 

m 

T = the temperature at the inside of the pipe wall. 

The ambient temperature is inputted with the materials properly card. 

The general portion of the NASTRAN program uses the load, stiffness 
and mass matrices to compute the displacement vector for each grid 
point. Force and stress recovery procedures are different for each type of 
element. In the general NASTRAN scheme, matrices for recovery of element 
forces are set up in one subroutine and actual force recovery takes place 
in another subroutine. In the first subroutine, the matrices required to 
determine forces at end (a) due to displacements of end (a) and 
displacements of end (b) are generated, along with matrices which include 
the original forces due to temperature expansion. These matrices are 
defined by the equations 
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(26) 


and 


[SJ 


,R 


(a) 


a-, = [K aal ^ T eb 3 ^ ^ 

b X b 


/R 


(b) 


fi = ^ f T eb 3 ^ ^ 

6x6 


(27) 


PS T» 

6x1 


C 


X K " X , 

b a 


y b - y a 

o ) 


0 
0 

v 0 


I 

J 


Once the force recovery matrices have been determined, the forces 
at end (a) can be calculated by the equation 

(P a ) = [S a ] <P A } ♦ [S b ] (V + <V x (T -l m ) + {P D } (29) 

where 

{ P } = { P , P , P , M. , M. , M }, force vector at end a, 
a a R a 0 a Z a R a e a Z 

T AMB = am bi ent temperature, 

{P D> ■ [K aal , {Wfl, 6 x 1 

D X D 

and 

{ DFA} = displacement components at end (a) corresponding to enforced 
deformation. 

The normal force is redefined so that tension is positive: 


F = - P 
r G r 0 ' 
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Once the forces at end (a) have been found. Equation (7) can be 
used to determine the forces at end (b) of the element. 

The stresses to be calculated for the curved pipe element are: 
For the element: 


Hoop stress: 


a H = 


c E P r d 

zt 


(30) 


Linear temperature gradient stress: 


C F E a (T 0 - T,) 

2 U~v) 


Nonlinear temperature gradient stress: 

E a A T 0 


(l -vT ; 


At each end: 


Axial stress: 


C A F 0 


°A = ~T 


Torsional shear stress: 


c b M e d 


Bending stress: 


c n 1 7 < M r) 2 + < M Z )2 


Lateral shear stress: 


t f = 


c c y'(p R ) 2 + (p z ) z 
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where 


^ b C, D, E, F = stress ^tensity factors input by the user, 

E = modulus of elasticity, 
a = coefficient of thermal expansion, 
d = outer diameter, 
t = wall thickness, 
v = Poisson's ratio, 

A = cross sectional area. 



< 
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Incorporation of Dummy Element Into NASTRAN 
The subroutines developed for a curved pipe become an operable 
reality in NASTRAN only after they have been coded in FORTRAN IV and 
only after several links in NASTRAN have been changed to allow for the 
use of dummy elements. A programmer familiar with NASTRAN data files 
and tables could eventually produce the necessary FORTRAN IV coding from 
scratch; however, a significant amount of time can be saved if the 
subroutines are patterned after an existing set of NASTRAN subroutines, 
approach also increases efficiency since it allows the programmer to 
take advantage of existing NASTRAN subroutines and accounts for NASTRAN 
idiosyncracies, such as the preference for recalculation of matrices 
rather than storage of the matrices. FORTRAN IV coding used for the 
curved pipe dummy element closely parallels the coding used in NASTRAN 
for the bar element. NASTRAN subroutines GMMATS (GMMATD for double 
precision) and TRANSS (TRANSD for double precision) were particularly 
convenient in making the necessary changes to bar element subroutines. 
The former, described in Articles 3.4.32 and 3.4.33 of Reference (8), 
is a general matrix multiplication and transpose subroutine, while the 
latter, described in Articles 3.4.37 and 3.4.38 of the same reference, 
is a subroutine which maps a vector from a local coordinate system to 
the basic coordinate system. 


4 


This 
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Communication between the subroutines and the functional module 


or other subroutines is simplified due to the similarity between the 
curved pipe element subroutines and the NASTRAN bar element subroutines. 

For all of the subroutines except the one which recovers element forces 
and stresses, it is necessary to retrieve information inputted on GRID, 

CDUM3, PDUM3, DEFORM, and MAT1 bulk data cards. Information on the 
DEFORM and MAT1 cards are stored and retrieved for the dummy element 
exactly as for the bar element. Information from the applicable GRID, 

CDUM3, and PDUM3 cards are stored in an orderly fashion in an array called 
the Element Connection and Property Table (ECPT). The format of the ECPT 
for all NASTRAN elements including dummy elements is described in Article 
2. 3. 8. 3 (and preceding articles) of Reference (8). The format of the ECPT 
for the dummy element depends on entries on the bulk data ADUMi card which 
defines the attributes of the dummy element. Generally, the format consists 
of: element identification number; all information on the element connection 

card after the property identification number; the material identification 
number; all information given on the element property card after the 
material identification number; the coordinate system identification 
number and basic coordinates for the grid points associated with the 
element; and the element temperature. Declarations in the bar subroutines 
were altered for the dummy element subroutines only to account for differences 
in the ECPT and in appropriate internal parameters. NASTRAN bookkeeping 
parameters were not changed. 
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The routine which recovers element forces and stresses (SDUM32) 
receives, as partial input, the output for the subroutine which sets 
up the element force recovery matrices. Most of the communication 
required between this subroutine and the functional module SDR2 
(Stress Data Recovery - Phase 2), is done with data block C0MM0N/SDR2X7/ 
described in Article 4.46.9 of Reference (8). The first 100 locations 
of this block are reserved for input parameters, the second 100 locations 
are for stress output parameters, and the third 100 locations are for force 
output parameters. The actual output format for dummy element forces and 
stresses is built into NASTRAN, rather than being set by the user. This 
is undesirable for the curved pipe element since the NASTRAN formats 
allow only nine forces and nine stresses per element as output, whereas 
a format allowing twelve forces and fourteen stresses is desirable in 
order to account for conditions at both ends of the element. The possibility 
of using WRITE statements in the force and stress recovery subroutine was 
considered and rejected because the dummy element output would then be 
located at random in the output of a NASTRAN job and would be difficult to 
interpret. The dummy formats could be changed using ALTER statements but 
this approach appears to be quite difficult. One alternative would be to 
output stresses instead of forces in the force output format since the pipe 
system analyst is more interested in stress output than force output. A 
second alternative would be to output forces and stresses at only one end 
of the element, since forces, and, therefore, stresses, can then be 
calculated for the opposite end. This is the approach presently being used. 
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Once the dummy element subroutines have been completed, a "linkedit" 
run is necessary to load them into the NASTRAN executable. The "linkedit" 
involves changes within NASTRAN itself. In this case the necessary 
"linkedit" programming was provided by the NASTRAN Systems Management Office, 
NASA-Langley Research Center. It was found that the "linkedit" used for 
one level of NASTRAN is not always completely acceptable for a different 
level . 

The most efficient approach to checking out a dummy element is to 
create "driver" programs which simulate NASTRAN with respect to each 
of the element subroutines so that trial input can be used to check 
each subroutine. All subroutines would be checked in this manner before 
the "linkedit" occurred. However, time limitations have prevented 
development of these "driver" programs for this effort, so present 
attempts at checking involve actual NASTRAN runs on a single element test 
problem. 

Once the dummy element subroutines have been debugged, application 
of NASTRAN with the dummy curved pipe element to the analysis of piping 
systems with curved sections can be compared to solutions from other sources. 
A NASTRAN formulation of Example #1 of Reference (9) has been prepared 
for comparison. Appendix II gives detailed information about the test 
problem and the test problem from Reference (9). 
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CONCLUDING REMARKS 


At this point significant advancement has been made in the 
development of a curved pipe element. The new element will simplify 
curved pipe analysis primarily because a single element is used 
rather than multiple elements, the geometry definition is simplified, 
and input/output are tailored specifically to pipe analysis. For 
generality, the flexibility factor and stress intensification factors 
will be inputted by the user. The NASTRAN user dummy element has 
proven to be a relatively simple means of instituting a new element, 
although format restrictions have caused problems. Use of subroutines 
from an existing NASTRAN element simplify communications with the 
NASTRAN functional modules but could produce errors in programming 
in sections where the two elements differ. The NASTRAN Systems 
Management Office provided the necessary linkedit routine and also 
provided assistance in removing inherent NASTRAN errors. Apparently 
the linkedit may require modification for different NASTRAN levels. 
Once the curved pipe element has been completely debugged, a straight 
pipe element should be considered so that piping analyses could be 
performed completely with appropriate NASTRAN input and output. 
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APPENDIX I 


UTILIZATION OF CURVED PIPE ELEMENT IN PIPING SYSTEM ANALYSIS 

The NASTRAN user can use the dummy curve pipe element almost as if it 
were a standard NASTRAN element. Some variations have been made from 
"standard" NASTRAN bulk data formats so that input will be as efficient 
as possible. Terns which may be dependent on the radius of curvature 
are restricted to the connection card, since the radius of curvature 
is likely to vary for different elements. Pin flags and offsets, 
conditions not expected to occur often in piping systems, were moved from 
connection cards to property cards to make room for the radius of 
curvature dependent terms. One ADUM3 card with non-varying entries must 
be included to describe the nature of the dummy element used. The required 
ADUM3, CDUM3, and PDUM3 cards are defined on the following pages. A sketch 
of the curved pipe element forces and geometry is given in Figure (4). 

Output for the element consists of forces and stresses at end (a) 
of the curved pipe (see Figure (4)). The forces given are, in order, 
the radial shear force (outward is positive), the axial force (tension 
positive), the out of plane shear force, and the moments in the element 
R, ip, and Z directions at end (a). Stresses given for end (a) are: 
hoop stress, linear temperature gradient stress, nonlinear temperature 
gradient stress, axial stress (tension is positive), shearing stress 
due to torsion, maximum bending stress, and shearing stress due to 
lateral loads. 
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The dummy element output formats provide headings for forces only as 
F1-F9 and headings for stresses only as S1-S9. Table (1) describes the 
force and stress output. Note that the outputted stresses have been 
pre-multi plied by stress intensification factors CA through CF, which 
were input by the user. 
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Figure (4). - Element Forces and Geometry. 
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BULK DATA DECK 


Input Data Card ADUM3 Dummy Element Attributes 

Description : Defines attributes of curved pipe element 

Format: 


ADUM3 

NG 

NC 

NP 

ND 


ADUM3 

2 

12 

14 

6 




Field : 

NG Number of grid points connected by DUM3 element (2) 

NC Number of additional entries on CDUM3 card (12) 

NP Number of additional entries on PDUM3 card (14) 

ND Number of displacement components at each grid point 

used in generation of differential stiffness matrix (6) 


Remarks : 

1. No differential stiffness matrix is actually generated. 

2. Entries on ADUM3 card for curved pipe elements will always be 

as shown above. 

3. Only one ADUM3 card should be present in the bulk data deck. 


i 
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BULK DATA DECK 


Input Data Card CDUM3 Curved Pipe Element Connection 

Description: Defines a curved pipe element using the dummy element 

capability 


Format and Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
L 

CDUM3 

EID 

PID" 

GA 

1 GB 

XI 

X2 

X3 

R 

abc 

CDUM3 

66 

21 

170 

! 180 



1 . 

4.5 

DC3 


+bc "1 

PI 1 

F 

CA 

CB 

CC 

CD 

CE 

CF 

1 

i 

i 

+C3 

I 10. 

3.2 

l 






1 


Field: 


EID 

PID 

GA, GB 
XI, X2 , 


R 


Unique element identification number (Integer > 0) 
Identification number of a PDUM3 property card (Integer > 0) 

Grid point identification numbers of connection points 
(Integer > 0; GA t GB) 

X3 Components of vector v at end (a), measured at end (a), 

parallel to the components of the displacement coordinate 
system for GA, to determine (with the vector from end (a) 
to end (b)) the orientation of the element coordinate 
system. The vector v lies in the plane of curvature of 
the element and is directed from end (a) to the inside 
of the curve and to the inside of the vector from end (a) 

to end (b). (Real, XI ^ + X2^ + X3 >0) 

Radius of curvature of the element (Real) 


PI Internal pressure (Real) 


F Flexibility factor (Real) 

CA CB Stress recovery coefficients for direct stresses, torsional 
CC* CD! stresses, shear stresses, bending stresses, hoop stresses, 

CF linear temperature gradient stresses, non-linear temperature 

gradient stresses (Real; default = 1) 


Remarks : 

1. Element identification numbers must be unique with respect to all 

other element identification numbers. 

2. See the sketch for a description of the curved pipe geometry. 
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BULK DATA DECK 


Input Data Card PDUM3 Curved Pipe Property 

Descri ption : Defines the properties of a curved pipe element 

~ " referenced by a CDUM3 card 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

PDUM3 

PID 

MID 


mum 

■MM 

mam 

warn 

TI ] 

abc 

PDUM3 | 

21 

10 



|| 






+bc 

PA 

PB 

M4I& 

e am 

wsmmm 

mm 

' Z3B ! 1 

+P3 







t i 

1 ! 


Fi el d : 


PID Property identification number (Integer > 0) 
MID Material identification number (Integer > 0) 
DO Outer diameter (Real) 


TH Wall thickness (Real) 

NSM Non-structural mass per unit (Real) 

length 

TO, TM, TI Outside, middle, and inside wall temperature (Real) 


PA, PB Pin flags for pipe ends (a) and (b), respectively, that are used 
to insure that the pipe cannot resist a force or moment 
corresponding to the pin flag at that respective end of the 
pipe. (Up to 5 of the unique digits 1-6 anywhere in the 
field with no imbedded blanks; integer > 0) (These degrees 
of freedom codes refer to the element forces and not global 
forces) . 

Z1A, Z2A, Z3A, Z1B, Z2B , Z3B Components of offset vectors w, and uj. , 

d D 

respectively, in displacement coordinate systems at grid 
points at ends (a) and (b) , respectively (Real or blank). 
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Remarks : 


1 . 

2 . 

3. 

4. 


PDUM3 cards my only reference MAT1 material cards since the 
curved pipe element is restricted to structural applications 


Note that separate property cards will be necessary when pin 
flags and/or offset are used. 


See the sketch for a description of the curved pipe geometry. 


If there are no pin flags or offsets, the continuation card may be 
omi tted. 
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TABLE (1) 


Curved Pipe Dummy Elements 
Element Printout Explanations 


Force Output 
Label 


EL ID 
FI 

F2 

F3 

F4 

F5 

F6 

Stress Output 


EL ID 

51 

52 

53 

54 

55 

56 

57 


Explanation 

[All forces at end (a)] 

Element identification number 

In plane shear force component (outward 
is positive) 

Axial force (tension is positive) 

Out of plane shear force component 
Bending moment about radial axis 
Torsional moment 

Bending moment about axis normal to plane 


[All stresses at end (a)] 

Element identification number 
Hoop stress 

Linear temperature gradient stress 
Nonlinear temperature gradient stress 
Axial stress (tension is positive) 

Shear stress due to torsion 

Maximum bending stress 

Maximum shear stress due to lateral loads 
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APPENDIX II 


TEST PROBLEMS 

The one element piping system shown in Figure (5) was used for 
initial tests of the dummy curved pipe element on NASTRAN. This problem 
was also solved by hand as a check on the NASTRAN run. Although the 
deflections and element forces found using the dummy element proved to 
be correct, the single point constraint forces were in error. The error 
can be traced to the transformation matrix involved in determining the 
global stiffness matrix from the element stiffness matrix. At this point 

the specific nature of the error is unknown. 

A second, more involved piping system problem was also coded for 
NASTRAN. This problem, shown in Figure 6, was taken from Reference (9) 
and provides a means of comparison of the NASTRAN solution with an 
available independent solution. Details of this problem are available 

in Reference (9). 
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R = 4.5" 


Outer diameter = 3.5" 

Modulus of elasticity = 

.229 x 10^ psi 

Moment of inertia = 3.017 in^ 
Shear Modulus = 1.15 x 10^ psi 
ijj = 90° 

Flexibility factor = 2 




Z 




./ i 




(basic coordinate system) 


X 


Figure (5). - Test Problem - Single Curved Element. 
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